Quantum Hall effects in layered disordered superconductors 
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Layered singlet paired superconductors with disorder and broken time reversal symmetry are 
studied. The phase diagram demonstrates charge-spin separation in transport. In terms of the 
average intergrain transmission and the interlayer tunnelling we find quantum Hall phases with spin 
Hall coefficients of a^ m = 0, 2 separated by a spin metal phase. We identify a spin metal-insulator 
localization exponent as well as a spin conductivity exponent of w 0.9. In presence of a Zeeman 
term an additional a^y 1 ™ = 1 phase appears. 
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The problem of quasiparticle transport and localiza- 
tion in disordered superconductors is of considerable in- 
terest in view of experimental activity on the high T c 
cuprates as well as theoretical realization that disordered 
superconductors provide new symmetry classes of ran- 
dom matrix theoryi Of particular interest is class C for 
which the Hamiltonian breaks time reversal symmetry 
but spin rotation invariance remains intact. Physically, 
it can be realized in materials consisting of singlet super- 
conductor grains in a magnetic field or else, by a super- 
conductor in the absence of magnetic field whose order 
parameter breaks time reversal invariance, such as d+id' . 
Class C can therefore be realized by high T c compounds 
where d wave pairing is well established. In fact, d + id' 
pairing has been suggested 2 , in particular in overdoped 
compounds or as field induced pairing 3 . 

Transport properties of random superconductors are 
unusual since a quasiparticle does not carry charge, be- 
ing screened by the condensate, while the singlet paired 
condensate does not transport spin. Furthermore, the 
gapless nature of d wave pairing with low lying quasi- 
particle excitations leads to a rich phase diagram in 2- 
dimensions (2D) with spin quantum Hall phases^ 5 -, spin 
insulators and spin metals&i; a metallic phase was also 
found for triplet pairing 8 . 

The usual quantum Hall system in two dimensions 
(2D), as well as its extension to 3-dimensional (3D) lay- 
ered system, have been studied by a network mode&ifl 
which consists of a lattice of nodes connected by links. In 
2D, the unidirectional propagation on links is described 
by random phases, corresponding to a group U(l), while 
transfer at nodes is controlled by a parameter which de- 
termines the critical point. The transfer matrix of the 
network model can be efficiently evaluated identifying 
the critical behavior, e.g. the localization exponent is 
vqh ~ 2.5. The 2D class C problem has recently been 
studied by a network model 4 where propagation on links 
of particle-hole spinors via the Bogoliubov-de Gennes 
Hamiltonian is realized by random SU(2) matrices. The 
quantized spin Hall conductance is shown to jump by two 
units at a critical point of a new universality class with 
a localization exponent vsd ~ 1.12; an exact mapping on 



a classical percolation problemi 1 * 1 ^ has found v%d = 4/3. 
The spin rotation invariance can be broken by having a 
different transmission for particles and holes^, e.g. a Zee- 
man term. The phase diagram has then 3 phases with 
quantum Hall values of 0, 1, 2, respectively and a local- 
ization exponent vqh of the usual U(l) theory. 

Experimental realization of this unusual spin transport 
depends on the ability to control deviations from the crit- 
ical point. In the usual quantum Hall effect this is con- 
trolled by the position of the Fermi energy relative to 
that of an extended state in a Landau band. In a su- 
perconductor the particle-hole symmetry fixes the Fermi 
energy at the middle of the gap, and the relative position 
of states is not directly tuned by the overall density. It 
was in fact suggested that changing the strength of dis- 
order can lead to quantum Hall transitions^, at least for 
weak breaking of time reversal symmetry. 

An important insight into the nature of d + id' super- 
conductors comes from studying their edge states^ 3 ^ 3 - 
which provide a realization of our network model and 
identify its parameters. In the d wave case a promi- 
nent zero bias anomaly 2 ^ has identified a surface state 
at zero energy. The d + id' case allows current carrying 
chiral states that split the zero bias anomaly as seen in 
the overdoped compounds^. The chirality of these edge 
states leads directly to quantized Hall conductanc e 5 ' 13 . 
In contrast, charge transport of superconducting grains 
is dominated by the randomness of the Josephson cou- 
pling between the grains; the phase correlation between 
grains is lost in 2D at a critical value of disorder, as shown 
in an XY model with random phase shifts^. 

In the present work we solve a network model for a 
layered 3-dimensional (3D) random superconductor. We 
find 3 phases which we identify as spin insulators with 
Hall coefficient and 2, respectively, and a spin metal 
phase. We also identify the localization exponent at the 
insulator-metal transition, v^d « 0.9. When the spin ro- 
tation symmetry is broken by a Zeeman term we find 
an additional phase with Hall coefficient of 1. The spin 
metal phase is a realization of a proposed phase with 
nonzero spin diffusion constant at zero temperature^. We 
identify the conductivity of the spin metal and find its 
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critical exponent to be v%d as well. Finally, we identify 
the physical parameter that controls criticality of spin 
transport, namely the average transmission of quasiparti- 
cles between grains. This then demonstrates spin-charge 
separation in the sense that their critical behavior rele- 
vant to transport is controlled by distinct parameters. 

The 3D network model consists of layers of 2D lattices. 
Each 2D lattice has nodes, which are connected by links 9 . 
Quasiparticles propagate unidirectionally on links, hence 
their transfer matrix is equivalent to the evolution op- 
erator, which for a singlet superconductor is an SU(2) 
matrix describing an Andreev process where particle and 
hole components mis^, 
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where <5i, 62, x are random, < Si, 62 < 2tt, < x < 1. 

The propagation between grains, i.e. at nodes of 
the network, is determined by the transmission between 
grains. At each node we have two incoming and two out- 
going links with particles and holes separately scattered. 
The transmission probability at a node is parameterized 
in the form Tq — [1 + exp(— 7re)] _1 , so that for e = the 
transmission equals reflection i.e. maximal mixing of all 
links. In the following we allow also a Zeeman parameter 
A which breaks spin rotation invariance. The transfer 
matrix across a node is then 



Vl + e -^±A/2) 

e -7r( e /2±A/4) 



,-7r(e/2±A/4) 



vT 



-7r(e±A/2) 



(2) 



where the ± sign corresponds to particle with spin-up or 
hole with spin-down, respectively. 

The final ingredient of the 3D network model are addi- 
tional nodesAfi connecting neighboring layers. Placing 2D 
layers one on top of the other allows for links belonging 
to neighboring layers to form nodes and quasiparticles 
scatter. The matrix describing this scattering is 
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Consider a system of size M x M x L where M is the 
number of links in one layer (with two channels per link) 
and L — > 00 is its length. For a given M the eigenval- 
ues of T^T, where T is the full transfer matrix, behave 
as e~ 2A " L , defining the Lyapunov exponents {A„}; the 
smallest positive one, A 1; defines the localization length 
£m = The M dependence of £,m/M identifies the 

phases: (i) a decreasing ratio corresponds to localized 
state, i.e. a spin insulator, (ii) a constant ratio corre- 
sponds to a critical state, and (iii) an increasing ratio 
corresponds to a spin metal. The phase diagram for 
the class C network model (with A = 0) is displayed in 
Fig. 1. Square boxes represent critical e cr (t) lines. The 
particle-hole symmetry of the superconductor ensures a 
degeneracy at the critical point e = t = 0, i.e. the Hall co- 
efficient changes by two units. Furthermore, in the clean 
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FIG. 1: Phase diagram of the 3D network model with spin- 
rotational symmetry. 



limit the Hall conductance has two unitsi 3 -; this corre- 
sponds to transmission To = 1, i.e. e > large. Hence 
there arc three distinct phases: Hall insulator with Hall 
conductance a xy — 0, spin metallic phase, and a quan- 
tized spin Hall phase with a xy = 2. 

The width Wit) in e of the metallic region increases 
with t, and is expected to behave asi^ W(t) ~ t 1 /"" 
where V2d is the localization length exponent in 2D. The 
argument is that for a 2D isolated layer the mean level 
spacing is ~ l/£fd w ith £2^ the 2D localization length. 
The states in each layer are concentrated along a perco- 
lation cluster of length £2^ and width of the edge state 
(a coherence length^) , i.e. normalized as ~ ^ad^'i * ne 
interlayer coupling is then ~ tf&d- At the mobility edge 
the mean level spacing is of the same order as the inter- 
layer coupling according to the Thouless criterion, hence 
e ~ t 1 l* J ' 2 < i . The curve in Fig. 1. represents the least 
square fit for the data, producing W ~ t 1 ' 1 ' 2 , which is in 
good agreement with our previous valued of V2d- 

The divergence of the localization length at e cr (£) iden- 
tifies the localization exponent of a spin insulator-metal 
transition in 3D. We have evaluated the critical expo- 
nent at the symmetric point t = l/v2 (larger t maps 
into a smaller t by rearranging layer indices) and found 
vzd = 0.91 on both the insulator and metal sides, as 
shown in Fig. 2; the same value is found for t = 0.1. 

Fig. 2 shows that in the metallic phase £m /M increases 
approximately linearly with M. It was proposed^ that 
this identifies the 3D conductivity as a xx ~ (e — £ C rY 3d - 
This derivation needs to be revised since in the 3D 
limit the conductivity involves many Lyapunov expo- 
nents A„ . The multichannel conductance is given byi&^£ 
g = £) n [l + cosh(A„L)] -1 . For a few channels M 2 < L 
the lowest Lyapunov dominates, but in the 3D limit 
M 2 ^> L the number of terms N e ff that contribute to 
g is large. In fact, for many channels the rigidity in the 
spectrum of T*T suppresses fluctuations in {A„} and one 
expects 17 A n = nAi. Hence g w N e ft sa l/(Ai-L). The 
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FIG. 2: Scaling of renormalized localization length for the 3D 
network model with spin-rotational symmetry at t — 1/V5 
showing an exponent v-^d = 0.91. The lower and upper 
branches correspond to the insulating e > e cr and metallic 
e < e cr phases, respectively. 

conductance has then the form 

g^N e ff^^-[a{e-e cr ) V3d M + b}. (4) 

This shows that the conductivity in 3D is indeed a xx = 
gL/M 2 ~ (e — e cr ) l/3d . On the critical line e = e cr the 
conductance is limited to the surface area and is ~ b. 

Consider next the A ^ case with broken spin rotation 
symmetry. At A = 2, e.g., the 2D system (t = 0) is 
critical at e cr — ±0.64 with a critical exponent vqh ~ 2.5 
of the usual quantum Hall system. At t ^ 0, we expect 
each of the critical states to split into two with a band of 
metallic states between them (as for A = 0), however it is 
not obvious wether the two internal critical curves merge 
or produce a new phase boundary. We find merging of 
these lines, producing a four-phase diagram as shown in 
Fig. 3. Both outer critical lines scale as t ~ |e — e C r\ VQH 
in agreement with the argument above. The inner curve 
is affected by both critical points and therefore deviates 
from this scaling form. 

Fig. 3 shows the existence of a new phase with a xy = 1 
which becomes metallic at very low values of t, e.g. 
t = 0.001 at e = 0. This feature can be traced to 
the rather large £m/M values for t = in the range 
—e cr < e < e cr . We note that at t — a single spin 
state becomes extended at e cr , while the other spin state 
becomes extended independently at — e cr . For t < 0.001 
these extended states produce two metallic bands which 
do not overlap, hence when the chemical potential is in 
between these bands the a xy — 1 phase emerges. When 
t > 0.001 these bands overlap and a a xy = 1 phase is not 
possible. We emphasize that there is a single metallic 
phase, hence in this phase the two extended spin states 
mix via interlayer coupling, unlike the situation at t = 0. 

We proceed to evaluate the localization length expo- 
nent for A = 2 and different t. For t = l/y/2 (maximal 



FIG. 3: Phase diagram of the 3D network model with A = 2 
showing an additional phase with a xy — 1. 

mixing of Eq. |2J we find v = 0.85 while for t = 0.1 we 
find v w 0.93, fairly close to v§&. It differs significantly 
from the value 1.45 found for the 3D U(l) systemic. This 
is consistent with our finding of a single metallic phase in 
which extended states with both spins are mixed. Thus 
at t = (extended states at ±e cr are decoupled) the sym- 
metry is reduced to U(l) producing the vqh exponent, 
while the t ^ metallic phase (mixed extended states) 
has the same symmetry for A = or A ^ 0. 

Based on these results we propose a 3-parameter scal- 
ing function near the multicritical point e = A = t = 

= f{e V2d M, A"M, t a M) (5) 

On the critical surface the scaling function is M inde- 
pendednt, hence for A = the critical line is e cr ~ t a / U2d . 
As discussed above, W(t) ~ e cr ~ t 1 > V2d , therefore a = 1. 

We demonstrate in Fig. 4 scaling along the line 
e = A = with a = 1. Note that when t is too large, 
approaching the symmetric point 1/V2, scaling is not 
expected. The analysis above yields a xx ~ t when ap- 
proaching the multicritical point. 

Finally we consider a realization of the parameter 
which drives the phase transitions. Correspondence with 
edge states of the d + id' system shows that the tran- 
sition is driven by the average value To of quasiparti- 
cle transmission between grains, which determines e via 
T = [1 + cxp(— 7re)] _1 . Following the suggestion that, at 
least for weakly broken time reversal symmetry, disorder 
may drive the transitions^, we have performed further 
simulations of the 2D network searching for the effect of 
randomness in e. We have found this randomness to have 
a negative scaling exponent, i.e. an irrelevant variable. 
In fact, the mapping to a percolation problem^ leads to a 
distribution in the percolation parameter which is indeed 
irrelevant. In our system where time reversal symmetry 
breaking is fully developed we do not expect disorder 
to drive a spin quantum Hall transition. In contrast, 
charge transport is dominated by the randomness of the 
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FIG. 4: Scaling of renormalized localization length for e 
A = as function of interlayer coupling. 



Josephson coupling between superconducting grains; the 
phase correlation between grains is lost in 2D at a critical 
value of disorder, as shown in an XY model with random 
phase shifts^. We conclude then that a superconductor- 
insulator transition for charge transport is disconnected 
from that of quasiparticle spin transport, realizing spin- 
charge separation in transport. 

As a concrete realization we consider two grains of d + 
id' superconductors with parallel edge states along an 
axis x. An impurity provides an intergrain scattering 
potential VaS(x), where a is a lattice constant. The right 
and left moving edge state ipn( x ) an d iPl(x) then satisfy 



ivd x ip R (x) + Va5(x)ip L (x) = Ei() R (x) 
ivd x ip L (x) + Va6(x)tp R (x) = Etp L (x) 



(6) 



velocity, v ss aA' for a (110) surface and v ss aA for a 
(100) surface^ where A and A' are the gaps of d x i_ y i 
and d xy , respectively, with A' « A. 

The transmission from an incoming ip R (x) to an out- 
going ^l(x) is readily evaluated as 



Tn = 



4(Va/2v) 2 
[1 + (Va/2v) 



(7) 



where E is an energy eigenvalue. Here v is the edge state 



Note that T has a maximum of 1 at Va/2v — 1 and 
decreases at large V [since then the matching of states 
near the impurity (E ±V) with the nearest levels on 
the edges (E ±v/a) is reduced]. The transmission 
needed for exhibiting an extended state, To — 1/2, is 
achieved at V « A' for (110) edges or V « A for (100) 
edges, i.e. a much weaker coupling in the former case. 

In conclusion, we have demonstrated spin-charge sepa- 
ration in transport: spin transport and related QH tran- 
sitions are controlled by the average intergrain transmis- 
sion while charge transport and superconducting corre- 
lation is controlled by the amount of disorder in the in- 
tergrain (Josephson) coupling. We show that interlayer 
coupling leads to a new spin metal phase and identify the 
localization exponent for the spin insulator-metal transi- 
tion as v-id w 0.9. The latter is also the spin conductiv- 
ity exponent when approaching the transition from the 
metallic side. 
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